Magnetization plateaus in the Ising limit of the multiple-spin exchange model on 

plaquette chain. 
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We consider the Ising spin system, which stems out from the corresponding Multiple-spin exchange 
(MSE) Hamiltonian, on the special one-dimensional lattice, diamond-plaquette chain. Using the 
technique e of transfer-matrix we obtain the exact expression for system free energy with the aid 
of which we obtain the magnetization function. Analyzing magnetization curves for varies values 
of temperature and couplin constants we found the magnetization plateaux at 1/3 and 2/3 of the 
full moment. The corresponding microscopic spin configurations are unknown by virtue of high 
frustration. 
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PACS numbers: 05.50.+q, 05.70.Fh 

The Heisenberg model is widely recognized as a lattice 
model for magnetism of materials. However, this model 
by no means is universal, because it based on several 
assumptions. One of these assumptions consists in the 
pair character of exchange interactions. This means that 
only the exchange processes of no more than two particles 
(nearest-neighbour or next-nearest-neighbour and so on) 
are taken into consideration: 
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Here P^- are the pair exchange operators, which imple- 
ment the transposition of two spin states in i-th and j-th 
sites of the lattice: 



(2) 



For the SU(2) spins and s = 1/2 the expression for P^ 
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<Ji are the Pauli matrices. The generalization of this 
picture is known since 60-s [Hand called the multiple- 
spin exchange (MSE) model. This model describes the 
magnetism of the system of almost localized fermions via 
the concept of many particle permutation. The general 
form of MSE Hamiltonian is 



H --V/ (-1Y P 
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Here P n denotes the n— particle cyclic permutation op- 
erator, J n is corresponding exchange energy and p is the 
parity of the permutation, which indicates how much 



pair transpositions contains in the given cyclic permu- 
tation. Many peculiar magnetic and thermodynamical 
properties of 3 He adsorbed on graphite surface can be 
understood only within the framework of MSE model 0] • 
Recently the significant role of MSE interaction was re- 
vealed in low-dimensional cuprate compounds [3|, what 
initiated the interest toward the two-leg spin ladders with 
four-spin cyclic interaction 0] ■ The MSE model by itself 
exhibits rich phase structure even at classical level 0]. 
Another interesting feature of MSE model is the pos- 
sibility of complex magnetic behavior, including such a 
phenomena as magnetization plateaux, which was estab- 
lished in the MSE model on triangular lattice [6j]. The 
study of the magnetization plateaux Q is the one of the 
main directions of nowaday investigations on macroscopic 
non-trivial quantum effects in condensed matter physics, 
which have a number of fundamental and applied impor- 
tance. Despite the purely quantum origin of this effect 
it was shown recently that, magnetization plat eaux can 
appear in the Ising spin systems also lij, llo|. exhibit- 
ing in some cases fully qualitative correspondence with 
its Heisenbegr counterpart [9(. The latter fact is very 
important because it can serve for more profound under- 
standing of magnetization plateaux physics, can provide 
it with new methods and can initiate a search of novel 
magnetic materials with huge axial anisotropy. 



I. THE MODEL 

We consider a one-dimensional lattice consists of cor- 
ner shared so-called diamond plaquettes (See Fig. [TJ. 
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FIG. 1: The diamond plaquette spin chain. 



The diamond plaquette is a square plaquette of 4 spins 
with nearest neighbor interaction and additional bound, 
connecting two opposite spins. If we consider the Hamil- 
tonian of MSE model and restrict ourself with the two-, 
three-, and four-spin exchanges, which is the simplest 
case, we arrive at he following Hamiltonian: 



(i, 3 ) p i 

h p = 2J °"i " °J + (°T ' °"2) (<T3 ' ^4) + (<Tl • CT 4 ) (<X 2 • <T 3 ) - (<Ti • <X 3 ) (<T 2 • <X 4 ) . 
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where J = J 3 — J 2 /2 and if = — J 4 /4. This model have 
been successfully applying for describing the properties 
of solid 3 He films. In our case we omit the off-diagonal 
part of the Hamiltonian[6]what is equivalent to the formal 



replacement of all spin operators by simple Ising variables 
which take values ±1. Putting the system, describing by 
the Hamiltonian[6]onto the lattice, depicted in Figure ??, 
we get the following model: 



- /3H = aiSiS i+ i + a 2 (s* + s i+ i) (U + r t ) + a 3 (Un + SiSi+iUn) + h (sj + U + 7?) , 

i 

(6) 



where 

01 = P ( J 3 - J 2 /2 - J 4 /4) , (7) 
« 2 =/3(J 3 /2- J 2 /2- J 4 /4), 

Here the Ising spins placed in the corners of diamond 
plaquette chain denoted by ti and r, and spins placed 
in the middle line by Si. The system allows the exact 



calculation of the partition function, which can be repre- 
sented a sa trace of the N-th power of the corresponding 
transfer-matrix [ill ] 

z = E e ~ fm = s p tN - ( 8 ) 

(s,t,r) 

The transfer matrix is 



2e ai+h (e 2a3 cosh(4a 2 + 2h) - 
2e~ ai (cosh(2/i) + 1) . 



2e~ ai (cosh(2/i) + 1) 
2e ai-h (- e 2a3 cos h(4a 2 - 2h) + e 2as 



(9) 



r 



As usually to calculate partition function we should ob- 
tain the eigenvalues of the transfer matrix, more precisely 



we need only the maximal eigenvalue, because in the ther- 
modynamical limit (N — » 00) only this one survives. So, 
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Zn = ^ lax The maximal eigenvalue is 



A = Acosh(h) + Bcosh(3h) 



. ^C fe cosh(2fc/i), 

\ k=o 



where the coefhcients are 



(10) 
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4e 2ai cosh(4a 2 ) 



8e 



-2ai 



C 2 = e 2a2 (2e~ 4ai + 2e 4a2 - e 4a3 ) 
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(ii) 

(12) 
(13) 

(14) 

(15) 
(16) 



Having the partition function we can obtain the free en- 
ergy of the system per one spin in the thermodynamical 
limit: 



Then , using the conventional thermodynamical relations 
rn = - (§L) , we obtain the magnetization per spin as s 
explicit function 



1 ,. logA^ 
- — hm 
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log A. 



(17) 



A sinh(/i) + 3B sinh(3/i) 



VELo co S h(2fe/i) 



3 ( A cosh(ft) + S cosh(3/i) + sj ELo C'fc cosh(2fc/i) J 



(18) 



Having this function we can draw the plots of the mag- 
netization processes for all finite temperatures and arbi- 
trary values of coupling constants. First of all, due to its 
geometry the system is highly frustrated in case of an- 
tiferromagnetic effective coupling constants. This means 
that even at T = the ground state of the system is 
disordered, because the arrangement of spins on lattice 
precludes satisfying all interaction simultaneously. How- 
ever, it is possible to choose such MSE coupling constants 
at which the frustration will be partially or entirely re- 



moved in the Ising limit. Among the variety of magneti- 
zation curves obtained for different sets of J2, J3 and J4 
the most remarkable is the one with two magnetization 
plateaux at m = 1/3 and m = 2/3 in the units of full 
moment(Fig. [J). At that region of coupling constants 
the system is highly frustrated and it is not so easy to 
determinate the microscopic spins configurations corre- 
sponding to these plateaux. Apparently they are some 
complex periodic structures with spatial period at lest 
3 for 1/3-plateau and 6 for 2/3-plateau. The investiga- 
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FIG. 2: The magnetization curve for Ising plaquette chain at T/J 3 = 0.17, J2/J3 = 1.5 and J4/J3 = 1.7. 



tion of the entire MSE model on diamond plaquette chain 
may change the picture obtained by us drastically. On 
the one hand, the appearance of other plateaux is possi- 
ble, on the other hand, the plateaux pertinent the Ising 
system might not survive in quantum case. The typical 
example is the simple S = 1/2 Heisenberg chain which is 



gapless, whereas the corresponding Ising system exhibits 
the plateau at m = 0. Let us mention also that appar- 
ently amon g th e overdoped RCuC>2 + x (R=Y, La, etc.) 
compounds 12| the ones are possible whose magnetic lat- 
tice is analogous to that considered here. 
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